On Hamiltonian flows whose orbits are straight lines 
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Abstract. We consider real analytic Hamiltonians on R"^ x K" whose flow depends linearly 
on time. Trivial examples are Hamiltonians H{q,p) that do not depend on the coordinate 
5 G R". By a theorem of Moser [5], every polynomial Hamiltonian of degree 3 reduces to 
such a (/-independent Hamiltonian via a linear symplectic change of variables. We show that 
such a reduction is impossible, in general, for polynomials of degree 4 or higher. But we give 
a condition that implies linear-symplectic conjugacy to another simple class of Hamiltonians. 
The condition is shown to hold for all nondegenerate Hamiltonians that are homogeneous of 
degree 4. 



1. Introduction and main results 

Polynomial Hamiltonians and maps have been studied extensively and for a variety of 
different reasons. Among other things, they constitute local normal forms for more general 
Hamiltonians and maps, and they provide a convenient testing ground for new ideas in 
dynamical systems. The restriction to polynomials also adds interesting algebraic aspects 
to the problem. This includes the possibility of classifying polynomial maps with a given 
property, and of decomposing them into simpler ones. 

The work presented here was motivated by a question, described below, that was 
raised in [8] concerning symplectic maps F with F — I a homogeneous polynomial. A 



differentiable map F : 



p2n 



— )■ 
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is said to be symplectic if 
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Here DF{x) denotes the derivative of F at x, and DF{xY denotes its transpose (as a 
matrix). If F — I is a homogeneous polynomial, then it can be shown [12] that F is the 
time-one map of a Hamiltonian H, and that H is affine-integrable as defined below. Thus, 
it is natural in this context to work with Hamiltonian systems. Moreover, some of our 
results do not require that the Hamiltonian be homogeneous, or even polynomial. 

Let if be a smooth function on M^"^. One of the basic facts from Hamiltonian me- 
chanics is that the vector field X = JVH defines a flow $ : (t, x) i— ;■ $*(x) whose time-t 
maps $* are symplectic. 

Definition 1.1. We say that a Hamiltonian H is afEne-integrable if its How $ is linear in 
time: 

$* = I-ftX, X = JVH, teR. (1.2) 

As we will see later, a Hamiltonian is affine-integrable if and only if the corresponding 
vector field X is constant along each orbit. That is, for all t G M we have 



X o{l + tX)=X . 



;i.3) 
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2 On Hamiltonian flows Avhose orbits are straight lines 

Polynomial maps F : R"^ — )■ M"^ with the property that X = F ~1 satisfies (1.3) are 
also called quasi-translations. They arise naturally in the study of singular Hessians [7]. 
To be more precise, the standard definition of a quasi-translation only requires that the 
identity (1.3) be satisfied for t = 1. However, this identity extends to t G N by induction, 
and further to t G C if X is a polynomial (since X o (l-\- tX) — X is a polynomial in t with 
infinitely many zeros). So I + tX is the time-t map for the vector field X. Differentiating 
the identity (1.3) with respect to t yields {DX)X = 0. Or equivalently, {X^V)'^£ = for aU 
linear functions i. This "local nilpotency" property is an alternative way of characterizing 
quasi-translations [10] and affine-integrable Hamiltonians [4,6]. 

In numerical analysis and physics, symplectic quasi-translations are also called jolt 
maps. They constitute the basic building blocks in the so-called Dragt-Finn factorization 
[2] of more general symplectic maps. This factorization has proved to be very useful in 
symplectic numerical schemes, including the simulation of Hamiltonian fiows in plasmas 
[4,9]. 

From a dynamical systems point of view, affine-integrable Hamiltonians are rather 
simple. Not only is the vector field X constant along each orbit, but its components Xj 
are Poisson-commuting invariants, as we will see later. So an affine-integrable Hamiltonian 
H is Liouville integrable, at least if it satisfies a suitable nondegeneracy condition. In 
addition, the geometry defined by the invariants Xj is quite restricted: 

Theorem 1.2. Let H be a real analytic affine-integrable Hamiltonian on M^". Then H 
and its vector field X are constant on the affine subspaces x + ra,nge{DX (x)) . If DX{x) 
has rank n then x + ra,nge{DX (x)) is a local level set for X . 

The only affine-integrable Hamiltonians that we have been able to find in the literature 
are all linear-symplectically conjugate to Hamiltonians of the form H{q,p) = K{p). The 
time-t map for such a Hamiltonian H is a shear: ^^(q^p) = [q + t'WK{p)^p^. 

Definition 1.3. We call H : R^"^ — )■ R a shear Hamiltonian if H is linear-symplectically 
conjugate to a Hamiltonian that does not depend on the variable q. In other words, there 
exists a linear symplectic change of variables U : R^"" — ?> R^"^ such that {HoU)(q,p) = K{j>) 
for some function K : R"" — )> R. 

Remark 1. As we will describe later, the linear map U in the above definition can be 
chosen both symplectic and orthogonal (as a matrix). 

One of our goals is to find an affine-integrable Hamiltonian that is not a shear, or 
to prove that there is no such Hamiltonian. Partial non-existence results can be obtained 
by restricting the class of Hamiltonians being considered. A trivial case: If n = 1 then 
line-orbits are necessarily parallel, so if iif : R^ — > R is affine-integrable, then there exists 
a rotation U of R^, such that [H o U){q,p) is independent of q. In other cases the Defini- 
tion 1.3 cannot be used directly. We shall give now give an alternative characterization of 
shear Hamiltonians. It applies to Hamiltonians that are regular in the following sense: 

Definition 1.4. We say that a Hamiltonian H is regular if there exists a point x where 
X{x) belongs to the range of DX{x). 
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Notice that homogeneous Hamiltonians of degree > 2 are regular, since their vector 
field vanishes at the origin. The following theorem was proved in [8] for Hamiltonians that 
are homogeneous polynomials of degree > 3. 

Theorem 1.5. A regular real analytic function H on M?'^ is a shear Hamiltonian if and 

only ifDX{x)DX{y) = holds for ah x,y E M^". 

To continue our discussion of special cases, assume that H is regular and affine- 
integrable. Then (DX)'^ = 0, as we will see later. In particular, if H is quadratic then 
Theorem 1.5 implies that H is a shear Hamiltonian. The cubic case is covered by a result 
of Moser [5] on quadratic symplectic maps on M^". It states that every such map F admits 
a decomposition F = AoS oL into three simple symplectic maps: an affine map A, a shear 
S{q,p) = {q + s{p),p), and a linear map L. It is not hard to show that this result implies - 
and is essentially equivalent to - the statement that every homogeneous affine-integrable 
Hamiltonian of degree 3 is a shear. A direct proof is given in Section 2. 

This raises the question [8] whether every homogeneous affine-integrable Hamiltonian 
on M^" is a shear. Locally, much more is true: It is well known that every smooth Hamil- 
tonian H is of the form H{q,p) = K{p) in some local symplectic chart, near any point 
where the vector field does not vanish. The local conjugacy (chart) is nonlinear in general. 
But if all orbits for H are straight lines with constant velocity, as is the case for affine- 
integrable Hamiltonians, and if X is constant on n-dimensional affine subspaces, then one 
might think that this conjugacy can be chosen to be linear. However, this is not true in 
general: 

Theorem 1.6. The following Hamiltonian on M.^ is affine-integrable but not a shear: 

H{q, p) = qipl + \/?>q2plpA + Pipl - V3p2P3P4 : g, p G R^ . (1.4) 

In addition, H is nondegenerate in the sense defined below. 

Definition 1.7. A real analytic vector field X on M^" is said to be nondegenerate if 
DX{x) has rank > n at some point x G M^". If X = JVH then we also say that H is 
nondegenerate. 

We would like to stress that this is the one and only notion of nondegeneracy that 
will be used in this paper. Notice that, if DX{x)'^ = 0, then the rank of DX{x) can be no 
larger than n. And if X is analytic, then the rank is constant outside some analytic set of 
codimension one. 

The example (1.4) belongs to a simple class of Hamiltonians that we shall now describe. 
Let < d < n. To simplify the description, we write q = (Q^q) and p = {P,p), where 
Q, P E M'^ and q,p E W^~'^. Consider a Hamiltonian of the form 

H{q,p)=K{p)+Q^V{p) + P^W{p), (1.5) 

with K : R*^-^ -^RandV^W : R'^''^ -^ R'^ differentiable. Notice that H does not depend 
on q, and thus p stays fixed under the fiow. Furthermore, the coordinates Q and P evolve 
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linearly (in time) under the flow. If q evolves linearly as well, then H is afline-integrable. 
As we will see later, this is the case if and only if 

wipyDVip) - vipyDWip) = o . (i.e) 

If d = 0, then P = Q = and p can be identifled with p. In this case, (1.5) becomes 
H{q,p) = K{p), so iif is a shear Hamiltonian. 

Remark 2. The Hamiltonian (1.5) can be written as the sum of Hi = Q^V{p) and 
H2 = K{p) + P^W{p). What makes this decomposition interesting is that the Poisson 
bracket {Hi,H2} = (V-ffi)^J(V-ff2) of Hi and H2 Poisson-commutes with both Hi and 
H2. 

Theorem 1.8. Let H be a nondegenerate real analytic afEne-integrable Hamiltonian on 
M?"^ . Then H is linear-symplectically conjugate to a Hamiltonian of the form (1 . 5) if and 
only ifDX{x)DX{y)DX{z) = for ah x.y.zE M^". 

Here, as in Theorem 1.2, a simple class of afline-integrable Hamiltonian is characterized 
by a nilpotency-type condition on the derivative of the vector field. This suggest there may 
be a natural hierarchy of such conditions, which characterize classes of increasingly complex 
afline-integrable Hamiltonians. 

We will prove that the condition on DX in Theorem 1.8 holds if if is a homogeneous 
polynomials of degree 4. As a result we obtain 

Theorem 1.9. Let H be a nondegenerate affine-integrable Hamiltonian on M^"^. If H is 
homogeneous of degree 4 then H is is linear-symplectically conjugate to a Hamiltonian of 
the form (1.5). 

As mentioned earlier, any symplectic map F : M^"^ — )■ M^"^ with the property that 
F — I is a homogeneous polynomial is the time-one map of an afline-integrable Hamiltonian 
[12]. Combining this result with Theorem 1.9 and using the decomposition described in 
Remark 2, we obtain the following factorization theorem. 

Theorem 1.10. Let F be a symplectic map on M^"^ such that F — 1 is homogeneous of 
degree 3 and nondegenerate. Then F admits a decomposition F = Fi o F4, where Fi and 
F4 are the time-one maps of two shear Hamiltonians. 

For the proofs of Theorems 1.2, 1.5, 1.6, 1.8, 1.9, and 1.10, we refer to Sections 3, 2, 6, 
4, 5, and 7, respectively. Some additional results can be found in Section 2. 



2. Basic properties 

In the remaining part of this paper we always assume that if is a real analytic Hamiltonian 
on M^". Furthermore, by a "homogeneous" Hamiltonian we always mean a homogeneous 
polynomial. 
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As is true in general, the Hamiltonian H is invariant under the flow that it generates, 
so {DH)X = 0. Assuming that H is affine-integrable, X o (I -\- tX) = X, and thus 
{DX)X = 0. Furthermore, $* = I + tX is symplectic, which by (1.1) yields 

J + t [{DXyj + 3DX] + t^ [{DXyjDX] = J . (2.1) 

Using that the terms of order t and t^ have to vanish separately, we get [DXy = 0. 
Differentiating the identity X o (l + tX) = X yields [DX o (I + tX)]{l + tDX) = DX. 
Multiplying on the right by (l — tDX) and using that {DXy = 0, we flnd in addition that 
DX o (I + tX) = DX. In summary, we have the following 

Lemma 2.1. Let H be an afEne-integrable Hamiltonian. Then the functions H and X 
and DX are constant along every orbit. Furthermore, {DH)X — and {DX)X = and 
[DXf = 0. 

An equivalent formulation of these facts can be given in terms of Poisson brackets. 
Assume that H is afflne-integrable. Let i.{x) = (Ju^x for some vector u G M^"". Then 
{£, H} is the directional derivative of H in the direction u, which we denote by duH. Being 
linear in the coordinate x, i evolves linearly in time, so {£, H} = duH is invariant under 
the flow. This implies the flrst identity in 

{duH,H} = Q, {dud^H,H} = 0, {duH,d^H} = 0. (2.2) 

The second and third identities are obtained from the flrst by applying a derivative 9„. 
This yields {dudyH, H} + {duH, d^H} = 0, and the two terms have to vanish separately 
since the flrst is symmetric in (m, v) and the second antisymmetric. If H is nondegenerate, 
then (2.2) shows that n of the vector fleld components Xj constitute a maximal set of 
Poisson-commuting invariants. So H is Liouville integrable, as mentioned earlier. 

Let X be a flxed but arbitrary point in M?"^. In the canonical splitting M?"^ = W^ x M"-, 
we can represent the derivative of X and the Hessian H[(x) = VDH{x) as 2 x 2 matrices 
whose entries are n x n matrices, 



DX{x) 



Z{xy A{x) 
-B{x) -Z{x) 



[X 



B{x) Z{x) 
Z{xy A{x) 



(2.3) 



Given that X = JVH, we have DX = JH. Since EI is symmetric, so are A and B. In 
the case of an afflne-integrable Hamiltonian, AZ is symmetric as well, as a result of the 
identity {DXf = 0. 

Lemma 2.2. A regular Hamiltonian is of the form H(q,p) = K(p) if and only if Z{x) = 
and B{x) = for all x. 

Proof. The necessity of the conditions Z = and B = is obvious. Assume now that they 
are satisfled. Let xq = (qo^Po) be a point where X{xo) = (Vp-ff(xo), —VqH^xo)) belongs 
to the range of DX{xo). At this point we have VqH{xo) = 0. Given that D^H = 
by assumption, this implies that the function q i— )■ H{q,po) is constant. Furthermore, 
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DpH{q,p) does not depend on g, since DqDpH = 0. So H{q,p) is independent of q as well, 
since H{q,p) = H{q,pQ) + f^ DpH{q,po + sv)vds with v = p — pq. qed 

Lemma 2.3. Let H be an afEne-integrable Hamiltonian. Given x G M?^, there exists an 
orthogonal symplectic 2n x 2n matrix U, and a diagonal n x n matrix A, such that 



U-^DX{x)U 



A 




u^mx)u 




A 



(2.4) 



Proof. Let d be the rank of M = E[(x). Let {ui,U2, ■ ■ ■ ,Ud) be an orthonormal set of 
eigenvectors for the nonzero eigenvalues of M. Since MJM = 0, the vectors Jwj are 
eigenvectors of M for the eigenvalue 0. Consider first the case d = n. Let U be the 2n x 2n 
matrix whose columns vectors are Jwi, . . . , J'Unj'Ui, . . . , -u^, in this order. Clearly, U is 
orthogonal and U^ MU diagonal. A simple computation shows that U is symplectic. 

If (i < n, consider the orthogonal projection P onto the span of J-ui, . . . , J-u^, tti, . . . , Ud- 
Then P commutes with both M and J. So we can choose an orthonormal set of vectors 
iud+i, • • • , Jwn, Ud+1, . . . ,Un in the null space of P and define U as above. qed 

The same construction can be used to give a 

Proof of Theorem 1.5. The necessity of the condition DX{x)DX{y) = is obvious. 
Assume now that this condition holds, for all x,y E M^"^, and that H is regular. 

First, we show that H is affine-integrable. By regularity, there exist xo,w G W^ such 
that DX{xq)w = X{xq). Thus X(x) = DX{xo)w + f^ DX{xo + sv)vds, for any given 
X G M^"^, where v = x — xq- This shows that DX(x)X(x) = for all x, which implies that 
H is affine-integrable. 

Let (-ui, ^2, • • • , Ud) be an orthonormal basis for the subspace spanned by all vectors 
M.{y)z with y,z E R^"^. Then EI(a;)J'Uj = for all x and all j. Defining U as in the proof 
of Lemma 2.3, we obtain (2.4) simultaneously for all x. (The matrix A can depend on x 
and need not be diagonal.) So {H o U){q,p) is independent of q by Lemma 2.2, implying 
that i7 is a shear Hamiltonian. qed 

As a corollary we obtain 

Theorem 2.4. [5] Every afEne-integrable Hamiltonian H that is homogeneous of degree 
3 is a shear. 

Proof. By (2.2) we have {d^H, H} = for /c < 2. The same holds for /c > 3 since H is 
of degree 3. It follows that {H{. + w), H} = for all u. Or equivalently, X{x)^JX{y) = 
for all X and y. From this we get DX{x)DX{y) = by differentiation, and the assertion 
follows from Theorem 1.5. qed 

The matrix U described in Lemma 2.3 is both symplectic and orthogonal. This means 
that U^JU = J and If^U = 1. As a result, we also have JU = UJ. In fact, any two of the 
three properties imply the third. This is know as the 2-out-of-3 property of the unitary 
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group U(n) = 0(?i) fl Sp(2n,M) fl GL(?i, C). The complex structure here is given by the 
matrix J, and the equation JU = UJ simply says that U is "complex" . Using the properties 
If^U = I and JU = UJ, any matrix U G U(n) can be written as 



U 



S T 
-T S 



S^S + T^T = l, S^T = T^S. 



(2.5) 



We will refer to such a 2n x 2n matrix as being unitary. The n x n submatrices S and T 
will be referred to as the real and imaginary parts of U, respectively. 

Concerning the claim in Remark 1, we note that any symplectic matrix M can be 
written as a product M = UAN, where U is unitary, A positive diagonal, and A^ unipotent 
upper-triangular. This is the standard Iwasawa decomposition [3]. If i7 is a Hamiltonian 
such that {H o M){q,p) is independent q, then (H o U){q,p) is independent of q as well. 

By Lemma 2.3, the Hessian W{x) of an affine-integrable Hamiltonian H is always of 
the form 



H(x) 



S T 
-T S 




A 









TAT^ TAS^ 
SAT^ SAS^ 



(2.6) 



for any given x G M^", where the matrix A can be chosen to be diagonal. This representa- 
tion is unique if M.{x) has n distinct nonzero eigenvalues, and if the diagonal elements of 
A are required to be in some prescribed order. 

If we do not require that the matrix A be diagonal, then we could replace S, T, and 
A in (2.6) by SV, TV and V~'^AV, respectively, where V can be any orthogonal n x n 
matrix. This fact is used in the lemma below. 

Example 3. Let M be an to x n matrix of rank m < n, and let / : M"^ — ?> M" be real 
analytic. Assuming Mf = 0, the equation q = f{Mq) defines a fiow on R"^ that is linear 
in time: $*((3'o) = Qo + tf{MqQ). This is similar to the flow considered in [11, Lemma 5]. 
It extends to a Hamiltonian flow on M^"^, as does every flow on M"^. The Hamiltonian is 
H{q,p) = p^ f{Mq). Using Theorem 1.5, is is easy to check that H is a, shear Hamiltonian. 
In fact, H can be triviahzed exphcitly: If we set T = M^{MM^)-^M and 5 = I - T, 
then (2.5) deflnes a matrix U G U(n), and we get {H o U){q,p) = H{p,p). 

Definition 2.5. We say that M.{x) is in semi-normal form if Z{x) = and B{x) = 0. 

Lemma 2.6. Let H be an affine-integrable Hamiltonian and x G M^"^. If A{x) is nonsin- 
gular then U'^M.{x)U is in semi-normal form for the matrix 



U = exp 



9 
-9 



9{x) = tan-i {({x)) , ({x) = Z{x)A{x)-^ . (2.7) 



Proof. Deflne ( = C{x) as above. A comparison with (2.6) shows that ( = TS~^. The 
conditions in (2.5) on 5" and T imply that ( is symmetric, and that 5*^(1 + C'^)S = I. 
Since SV = \S\ for some orthogonal matrix V, we can choose S" to be a positive definite 
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symmetric matrix. The choice is then unique: 5" = (I + C^)~^ • Setting 9 = tan~^(C) we 
obtain 5" = cos(6') and T = (^S = sin(6'), which leads to the expression (2.7) for U. qed 

This offers another way of checking whether H is a, shear Hamiltonian. First, we 
note that a nondegenerate affine-integrable Hamiltonian H is regular: X{x) belongs to 
the range of DX{x) at every point x where DX{x) has rank n, since DX{x)X{x) = by 
Lemma 2.1. 

Lemma 2.7. Let H be an affine-integrable Hamiltonian. Assume that A{xo) is nonsingular 
at some point Xq . Then H is a shear if and only if C, is constant near xq . 

Proof. First, assume that C, is constant near xq. So near xq, the matrix U in Lemma 2.6 
is independent of x, and H oU is iiv semi- normal form for a fixed unitary matrix U . By 
analyticity, this property extends to all x G M?^. Furthermore, ]H[(xo) has rank n, as (2.6) 
shows, implying that H is regular. So iif is a shear Hamiltonian by Lemma 2.2. 

Conversely, assume that {H o U){q,p) is independent of q for some linear symplectic 
matrix U. Then U can in fact be chosen unitary, as was shown the proof of Theorem 1.5. 
If 5" and T are the real and imaginary parts of U, as defined by (2.5), then we have 
Z{x)A~^{x) = TS~^ at every point x where A{x) is nonsingular. qed 

For completeness, let us mention that there is an alternative representation of M{x) 
via the shear map {q,p) i-> (q^p + Cl)- Assuming that H is affine-integrable and A = A{x) 
nonsingular, 

"I cl [0 oi ri oi _ \cAc cA 

1 ^ c I ~ ^C A 



m{x) 



(2.8) 



This shear is not unitary. But it is symplectic, since ( is symmetric. Furthermore, these 
shear maps form a group. Notice also that B{x) = Z{x)C{x). So the condition B{x) = 
in Definition 2.5 is redundant if H is affine-integrable and A{x) nonsingular. 



3. Invariant afRne subspaces 

In this section we give a proof of Theorem 1.2 and some related results. It is always 
assumed that H is affine-integrable and real analytic. 

Besides the flow $ for the Hamiltonian H, consider also the flows \E'j for the Hamilto- 
nians djH, where djH denotes the j-th partial derivative of H. By standard ODE results, 
\1/*- (x) is well deflned for all times t in some open neighborhood of zero in C (which may 
depend on x). By (2.2) the flows \l/j commute with each other and with $. So the flow 
\1/^ for dwH is given by 

Again, \E'^(a:) is well deflned for all w in some open ball B{x) C C^"^ centered at the origin. 
Furthermore, the "group property" \l/"(\l/^(x)) = \l/""'""'(x) holds whenever w, w+u G B(x) 
and ue S(^"'(x)). 

Lemma 3.1. Let u G M?^ . If the derivative of duH vanishes at some point x, then it 
vanishes at "^"^{x) for every w G B(x). 
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Proof. Define {G)t = G o ^i"" for any function G on R^^. Let now G = duH. Tfien 

Q = {d,{G,d^H})^= {{d,G,d^H} + {G,d,d^H})^ 
= {{djG)u d^H} + {G, {djd^H)t} . 

Here, we fiave used tfiat dwH and G are invariant under tlie flow \1/"', and tliat tlie maps 
\1/*™ are symplectic. Tfius, we liave 

j^{djG)t = {{djG)t, d^H} = -{G, {djd^H)t} = - J^id^G) J ^,rdr {djd^H)t . (3.3) 

(J,T 

Due to tlie factors d^^G that all vanish at x, the value {djG)t{x) is independent of t and 
thus (9jG)(\l/*™(x)) = {djG)t{x) = {djG){x) = 0, for all t in some open neighborhood of 
zero. The assertion now follows from the above-mentioned group property of \1/ and the 
analyticity of G. qed 

Corollary 3.2. Let u G M". If Hw vanishes at some point x, then M.u vanishes at \E'™(x) 
for every w G B{x). In other words, the nuU space (and thus the range) of HI is invariant 
under \E'^. 

Notice that the same holds for DX = JH. 

Corollary 3.3. Let x G M?"^ and R{x) = rsinge{DX{x)). Then \E'™(x) belongs to the afRne 
space X + R{x) for all w G B{x). Furthermore, w i— )> \E'™(x) is locally (near zero) invertible 
as a map from ,']R{x) to x + R{x). 

Proof. Consider the curve u{t) = \E'*™(a;)— x. Clearly u{0) belongs to R{x). The derivative 
u'{t) = DX{'^^^(x))w belongs to the range of DX{'^^'^{x)), which agrees with R{x) by 
Corollary 3.2. Thus, u{t) belongs to R{x) whenever tw G B{x). Since M.{x) is symmetric, 
DX{x) is invertible as a map from JR{x) to R{x). Thus, by the implicit function theorem, 
the same holds locally (near zero) for the map w i— ?■ \E'"'(x) — x, whose derivative at w = 
is DX{x). QED 

Proof of Theorem 1.2. For each y in x + R{x) there exists an open neighborhood 
By oi y in X + R{x) that is included in the range of fy : w ^-^ \E'™(y). This follows from 
Corollary 3.3. The vector field X is constant on each By since each component Xj is 
invariant under the flow \E'^. Similarly for H. Furthermore, the open sets By cover the 
afiine space x + R{x), and since this space is connected, it follows that X and H are 
constant on x + R{x). 

Assume now that DX{x) has rank n. If y = x + u + v, with u G R{x) and v G R{x)-^, 
then 

X{y) = X{x) + DX{x + u)v + 0{\vf) . (3.4) 

If y 7^ x is sufficiently close to x then \DX[x -\- u)v\ is bounded from below by a positive 
constant times l^;!, so we have X{y) = X{x) if and only ify — x = u& R{x). qed 
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4. Proof of Theorem 1.8 

We will write the given nilpotency condition on DX in the form 



(a:")JH(x)JH(x') =0, 



x",x,x' em.^'' 



(4.1) 



It is straightforward to check that this condition is necessary for H to be linear-symplectic- 
ally conjugate to a Hamiltonian of the form (1.5). 

Assume now that if is a nondegenerate (and thus regular) real analytic affine-integra- 
ble Hamiltonian that satisfies (4.1). Consider a point xq = {qo,Po) where ]HI(xo) has rank 
n. By performing a unitary change of variables, if necessary, we may assume that H(xo) 
is in semi-normal form, 

^0 
A{xo) 



H(a:o) 



In the case where x' — x" — xq, the property (4.1) implies that D^H — 0, so that 



(4.2) 



[X 



Z{x) 
ZixY A{x) 



X G 



X,2n 



(4.3) 



Thus, our Hamiltonian H has to be of the form 



if(Q,p) = /C(p)+Q^V(p) 



(4.4) 



with /C : W^ — 7> M and V : W^ — ?> W^ real analytic. Another consequence of (4.1) is that the 
range of 

~ z{xyz{x'y z{xyA{^x') - a{x)z{x') ' 

Z{x)Z{x') 



je(x)JH(x') 



(4.5) 



is contained in the null space of M.{x"), for every x" . If we take x" — xq then this implies 
that Z{x)Z{x') vanishes. But Z{q,p) = DV{p) and thus 



DV{p)DV{p') = , 



P,P e 



(4.6) 



Finally, multiplying (4.5) on the left by M.{x") and using that the result has to be the zero 
matrix, we find that Z(x"y A{x)Z{x') vanishes. In particular. 



DV{p"yK{p)DV{p') = 0, 



p\p,p' eW 



(4.7) 



where 



VDK, is the Hessian of /C. 



Let TZ be the linear span of all vectors V{p), with p G W^. Or equivalently, TZ is the 
linear span of all vectors DV{p)u with p,u & M". Here we have used that V{po) = 0. Then 
(4.6) and (4.7) imply that 



v{p + v) = v{p), ic{p + v) = ic{p) + Dic{py , ven 



{U 
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Let d be the dimension of 7^. If d = then there is nothing left to prove. Consider now 
the case where d > 0. Since V{p) = for all p G 7^, we also have d < n. Let TZ^ be the 
orthogonal complement of 7?. in M"^. Notice that q^V{p) vanishes whenever q G 7^^. 

Next we apply a unitary change of variables. To simplify the description, we rename 
current quantities by adding a subscript "old". The change of variables is (QomiPom) = 
(Sq^Sp), with 5" orthogonal, such that TZ = S~^7laM is the span of all vectors P = 
(pi, . . . ,pd, 0, . . . , 0). Then TZ^ is the span of all vectors p = (0, . . . , 0,pd+i, • • • ,Pn)- And 
V = S~^VomS takes values in TZ. From (4.8) we see that V{p) does not depend on P, and 
that }C{p) is an affine function of P. Setting Q = (gi, . . . , g^, 0, . . . , 0), the new Hamiltonian 
H is of the form (1.5). This concludes the proof of Theorem 1.8. 



5. Coramutators and quadratic functions 

The main goal in this section is to give a proof of Theorem 1.9 on quartic Hamiltonians. 
But some of the observations and computations apply to other Hamiltonians as well. 



The Hessian Il{u) at any point u G M^" defines a quadratic function x i— ?> 



l^T 



X M{U)X. 



The flow generated by this function is linear in time, since DX{u)'^ = 0. It is useful to know 
how these flows for different vectors u are related, as the hypotheses in Theorem 1.5 and 
Theorem 1.8 show. What simplifies the situation for homogeneous quartic Hamiltonians 
is that x^M.{u)x = u^M.{x)u = d^H{x), and if H is affine-integrable, then d^H commutes 
with H by (2.2). This fact will be exploited below. 

First we note that the Poisson bracket of two homogeneous quadratic functions 

F{x) = ^x^J-'x J G{x) = ^x^Qx, (5.1) 

is again a homogeneous quadratic function, 

{F,G}ix) = ^x^£x, J^= i(JJ^)(J^)-i(J^)(JJ^). (5.2) 

Here, J-", Q, and S are symmetric 2n x 2n matrices. The corresponding matrices JJ-", JQ, 
and JS belong to sp(2n,M). 

We also need to compute some double commutators, and not all functions involved 
are quadratic. To simplify the expressions, we use the operator notation |-f/^]-F = {F, H}. 
A straightforward computation shows that 

dud4Hj = {d^Hjdu + {duHjd^ + iHjdud^ + {dud^Hj 

Let now F be is a polynomial of degree < 2. Then d^dyF is constant and thus commutes 
with every function. If in addition F commutes with H, then (5.3) yields 

-Idud^HjF = IdyHjduF + Id^HjdyF , 
= d4dyH}F+ld^Hld,F. 

So far we have not used any properties of H other than differentiability. 
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Proof of Theorem 1.9. Assume now that H is homogeneous of degree 4. Let F be a 
polynomial of degree < 2 that commutes with H. Then the second identity in (5.4) implies 
that 

{d^'Hjdu'ld^HjduF = -ld,,HJldu'Hjd,duF = 0. (5.5) 

Assume in addition that H is affine-integrable. Then H satisfies (2.2). By Jacobi's identity 
for the Poisson bracket, if F and G commute with H, then so does |G]-F. In particular, 
{dudyHjF commutes with H. Using (5.4) and (5.5), we find that 

Idu'd^^HJldud^HjF = {Id^'Hjdu' + Idu'Hjd^^) {{d^Hjdu + lduH]d^)F = . (5.6) 

Since H is homogeneous of degree 4, the second derivatives of H are homogeneous quadratic 
polynomials. They commute with H by (2.2). So as a special case of (5.6) we have 

ldlH]ld^,HjdlH = 0. (5.7) 

Applying (5.2) with F = d^H and G = d^H yields 38 = 23M{u)3M{v) - 2JH(i;)JH('u). 
According to (5.7) we have JSJM{u) — JM{u)JS = 0, which implies that 

M{u)3M{v)JM{u) = , 'U,z;GM^". (5.8) 

Here we have used that EI('u)JHI(w) = 0. The identity (5.8) also be written as 

D^X{u,u,wym{v)D^X{u,u,z)=0, u.v.w.zeR'^''. (5.9) 

By Lemma 5.1 below, the same holds if each of the arguments u is replaced by a different 
vector in M^". In particular, we have (4.1) and thus DX {x")DX (x) DX {x') = for all 
x",x,x' G M^". The assertion now follows from Theorem 1.8. qed 

Lemma 5.1. Let X, Y, Z be vector spaces. Let (...): X^ -^ Y be a symmetric cubic form 
and : y^ — 7- Z be a symmetric quadratic form. Assume that (p, p, u) (p, p^v) =0 for 
all PjU,v E X . Then (wi, U2,us) (fi, ^2, vs) = for any Ui, Vj G X . 

Proof. Under the given assumption we have {p,p,p) {p,P,p) = 0. Now "differentiate" 
this identity twice: Replace p hj p + u + v, expand, and then collect all terms that are 
bilinear in {u,v). The result is 

12{p,u,v) {p,p,p) + 18{p,p,u) {p,p,v) = 0. (5.10) 

By assumption, the second term vanishes, so {p,u,v) {p,p,p) = 0. Differentiating this 
identity once we get 

(w, -u, v) (p, p, p) + 3(p, -u, v) (p, p,w) = . (5.11) 

Similarly, differentiating {p,p,u) {p,p,v) = once yields 

2{p,w,u) {p,p,v) + 2{p,p,u) {p,w,v) = 0. (5.12) 

Now (5.11) can be used to rewrite (5.12) as 

2 2 

--{v,w,u)Q {p,p,p) - -{p,p,p) {u,w,v) =0. (5.13) 

Or simplified, {p,p,p) (w, f, w) = 0. The assertion now follows by polarization. qed 
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6. The example from Theorem 1.6 

Consider the Hamiltonian (1.5), with the variables (Q, q; P,p) renamed to a; = (g, y;p, z), 

H{x) = K{z) + q'^V{z) + p'^W{z) , (6.1) 

where q^p eM.'^ and y^z E M""*^. The corresponding vector field X = [q^y-^p^z] is given 
by i = and 

q^W{z), p=-V{z), y,=d,K{z)+q^[d,V{z)\+p^[d,W{z)\. (6.2) 

Since i = 0, the Hamiltonian H is affine-integrable if and only if 

y,=q-'[d,V{z)]+p-[d,W{z)\ 

= w{zy[d,v{z)\-v{zy[d,w{-.^' 



is equal to zero for all j. This is precisely the condition (1.6). According to Theorem 1.5, 
if is a shear if and only if DX{x)DX{x') vanishes for all x and x' . Or equivalently, if and 
only if 

-x{xy:ix{x') = w{zyv{z') - v{zyw{z') (6.4) 

vanishes for all x and x' . 

It should be noted that the case (i = 1 is trivial: If v{s) = Vi{x + s{x' — x)) and 
w{s) = Wi{x + s{x' — x)) satisfy wv' — vw' = 0, then by the quotient rule of differentiation, 
the functions v and w are constant multiples of each other. So (6.4) follows from (6.3). In 
this case, H{q,p) can be made either independent of q via a change of variables [q^p] i— )■ 
[q^p + cq), or independent of p via a change of variables {q,p) ^-)■ (q + cp,p). Thus, if is a 
shear Hamiltonian ii d = 1. 

Notice also that, if the right hand side of (6.3) is equal to zero, then it remains zero if 
dj is replaced by 9? Thus, the right hand side of (6.4) is of the order \z — z'\^. This has 
motivated our choice of V and W below. 

The Hamiltonian (1.4) can be written as 

H{q,y;p,z) = qiVi{z) + q2V2{z) + piWi{z) +^2^^2(2) , (6.5) 

where g, y,p, z E M? and 

Vi{z) = z\ , Wx[z) = zl , V2{z) = V3zlz2 , W2{z) = -VSzizl . 

Let us compute the right hand side of (6.4), with z' replaced by w in order to simplify 
notation. If W2 = Z2 then 

[Wi{z)Vi{w) + W2{z)V2{w)] - [Vi{z)Wi{w) + V2{z)W2{w)] 

= Z2wl — ?>ZiZ2w\z2 — z\z2 + 3zfz2WiZ2 



3 (6-6) 

{wi-zifzl 



Wi — ?)ZiWi — Zl + 3ZiWi)Z2 
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This is clearly nonzero at some points, so H cannot be a shear. On the other hand, 

[Wi{z)d,Vi{z) + W2{z)djV2{z)] - [Vi{z)d,Wi{z) +V2{z)djW2{z)] =0 (6.7) 

holds for j = 1, due to the factor (t^i — zi)^ in (6.6). By symmetry, we have an expression 
analogous to (6.6) if wi = zi, with the cubic factor being {w2 — Z2)^ . So (6.7) holds for 
J = 2 as well. Thus, H is affine-integrable, as claimed in Theorem 1.6. 

A straightforward computation shows that the Hamiltonian (6.5) is nondegenerate: 
The Hessian E[(a;) has rank n = 4 whenever Z\Z2 7^ 0. Another noteworthy fact is that the 
matrix C,{x) defined in (2.7) depends on x only via the ratio z\l Z2- But the dependence is 
nontrivial, so by Lemma 2.7, this shows again that H cannot be a shear Hamiltonian. 

7. Elementary factorization 

A classical theorem by Jung [1] asserts that the group (under composition) of polynomial 
automorphisms of the plane R^ is generated by affine automorphisms and elementary shears 
(g,!?) I— )■ (q' + sd?),!?). No general result of this type is known in dimensions higher than 2. 
Theorem 1.10 covers the special case of symplectic maps F = I + X, with X homogeneous 
of degree 3. Its proof is based on the following observation. 

Lemma 7.1. Liet Hq be a polynomial afEne-integrable Hamiltonian of the from (1.5). 
Write Ho = Hi+ H2, where Hi = Q^V{p) and H2 = K{p) + P^W{p). Then Hi and H2 
Poisson-commute with if 3 = ^{Hi^ H2}. Furthermore, Hi and H4 = H2 — Hs are shear 
Hamiltonians, and the corresponding time-one maps satisfy 

$]^„ = $]^^ o $1,^ . (7.1) 

Proof. A straightforward calculation yields H^ = ^V{p)^W{p). So each of the Hamiltoni- 
ans Hj is of the form (1.5). Here, and in what follows, < j < 4. In addition, Hj satisfies 
the affine-integrability condition (1.6). Thus, the adjoint map {Hj} : G 1— )> {G, Hj} has the 
following nilpotency property: If / is any polynomial, then {HjJ / = for sufficiently large 
k. Furthermore, Hj commutes with if 3, since if 3 only depends on the variable p, while 
Hj is independent of the variable q. Thus, by the Baker-Campbell-Hausdorff formula, 

gt[Ho] J ^ gt[HiIgt[H2]g-t2[H3l J ^ gt[Hi]gt[H2-tH3] J ^ ^7_2) 

for every polynomial /. To be more precise, (7.2) is an identity for formal power series. But 
due to the above-mentioned nilpotency property, only finitely many terms of the series are 
nonzero. So (7.2) holds as an identity between polynomials. Using that / o $^ = e^^-^^^f 
for any polynomial /, we obtain (7.1) from (7.2). 

Let j > 1. Then the vector field Xj = JVHj satisfies Xj{x)^JXj{x') = for all x 
and all x' , as can be seen from (6.4). This shows that Hj is a shear Hamiltonian. qed 

Remark 4. The time-one map for if4 is an elementary shear, $^ {q, p) = (Q+V/i4(p), p), 
where h4{p) — H^^ip^p). The time-one map for ifi is unitarily conjugate to an elementary 
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shear Si{q,p) = {q + Vhi{p),p). A straightforward computation, similar to the one in 
Example 3, shows that hi{p) = Hi{p,p). 

Proof of Theorem 1.10. Let F be a symplectic map on M^"^ such that X = F — 1 is a, 
homogeneous polynomial of degree m >2. First, we prove that X — 3VH for some affine- 
integrable Hamiltonian H, using the same arguments as in [8] . The symplecticity condition 
(1.1) implies that X satisfies the equation (2.1) for t = 1. In this equation, the terms in 
square brackets have to vanish separately, since they have different degrees of homogeneity. 
The first of the resulting identities implies that the derivative of JX is a symmetric matrix. 
Thus, by the Poincare Lemma, JX is the gradient of a function —H. The second identity 
implies that (DX)'^ = 0. Thus, {DX)X = 0, since X{x) = m~^DX{x)x by homogeneity. 
This shows that the fiow for X is linear in time. In conclusion, F is the time-one map of 
an affine-integrable Hamiltonian H that is homogeneous of degree m-\-l. 

Consider now to = 3, and assume that X is nondegenerate. By Theorem 1.9, there 
exists a linear symplectic map U on M^", such that Hq = H o U~^ is a Hamiltonian of the 
form (1.5). In fact, U is unitary, as seen in Section 4. Using Lemma 7.1, we have 

F = ^\j = U-^o ^\j^^ oU = U-'o $1,^ o^\,^oU = ^jj^^u ° ^kou , (7.3) 

where Hi and H4 (and thus Hi o U and H^ o U) are shear Hamiltonians. qed 
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